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A general Hamiltonian theory for the adiabatic motion of relativistic charged particles confined 
by slowly-varying background electromagnetic fields is presented based on a unified Lie-transform 
perturbation analysis in extended phase space (which includes energy and time as independent 
coordinates) for all three adiabatic invariants. First, the guiding-center equations of motion for a 
relativistic particle are derived from the particle Lagrangian. Covariant aspects of the resulting 
relativistic guiding-center equations of motion are discussed and contrasted with previous works. 
Next, the second and third invariants for the bounce motion and drift motion, respectively, are 
obtained by successively removing the bounce phase and the drift phase from the guiding-center 
Lagrangian. First-order corrections to the second and third adiabatic invariants for a relativistic 
particle are derived. These results simplify and generalize previous works to all three adiabatic 
motions of relativistic magnetically-trapped particles. 

PACS numbers: 45.20.Jj, 52.27.Ny, 94.05.-a 

I. INTRODUCTION 

The concept of the adiabatic motion of a charged particle in magnetic fields is important to research in space 
plasma physics and fusion physics l|, |2|. Depending on the confining magnetic geometry, a particle may display 
three quasi-periodic or periodic motions. The fastest of these three motions is the gyromotion about a magnetic 
field line (with frequency w g ). The second motion exists when a particle bounces along a magnetic field line between 
two mirror points (with frequency Wb), because of nonuniformity along magnetic field lines. The slowest motion is 
the drift motion across magnetic field lines (with frequency u>d) caused by perpendicular magnetic gradient-curvature 
drifts. In space physics (and especially in radiation-belt physics [3(), these frequencies are widely separated such that 
uj g : LOh ■ oJd ~ e ■ 1 : e> where e <C 1 is a small dimensionless ordering parameter to be defined below. Associated 
with each periodic orbital motion, there exists a corresponding adiabatic invariant. We use /x, Jb and Jd for the three 



invariants constructed in this work; to be consistent, we may also use J g = {mc/q)n, where m is the particle's rest 
mass and q its charge. 

The theory of the adiabatic motion of charged particles in electromagnetic fields has been well developed by Northrop 
[jj. However, the non-Hamiltonian method used by Northrop resulted in dynamical equations that do not possess 
important conservation properties like energy conservation in static fields, because of the absence of higher-order terms 

n 

from Northrop's equations. In later work, Littlcjohn 14] used a noncanonical phase-space transformation method, 
based on Lie-transform perturbation analysis, to obtain the Hamiltonian formulation of guiding-center dynamics for 
nonrelativistic particles. By asymptotically removing the dependence on the gyrophase, the first invariant J g = 
(mc/q) /j, is obtained from the guiding-center Lagrangian by Noether's theorem. The resulting Hamiltonian guiding- 
center equations of motion conserve total energy for motion in static fields. In the present work, we use the Lie- 
transform perturbation analysis to develop a systematic Hamiltonian theory for relativistic guiding-center motion 
in weakly time-dependent electromagnetic fields. Our relativistic guiding-center equations of motion are expressed 

n n 

in semi-covariant form [5J, which simplifies the previous work by Grebogi and Littlejohn 6( (who extended their 
relativistic guiding-center equations to include ponderomotive effects associated with the presence of high-frequency 
electromagnetic waves) and generalizes earlier work by Brizard and Chan 7[ (who considered guiding-center motion 
of a relativistic particle in static magnetic fields). 

The derivation of relativistic guiding-center dynamics begins with the removal of the gyrophase dependence from 
the particle phase-space Lagrangian. Since the condition for these periodic motions to exist is that the time variations 
of the forces a particle experiences should be slow compared to the particle's motion, we assume first that the 
electromagnetic fields vary on the drift timescale. Thus we shall construct the first and second adiabatic invariants 
from the particle's motion. While this ordering is not the most general case, it is the most common one in practice 4|. 
This procedure gives us the reduced six-dimensional guiding-center Lagrangian and the first invariant J g . Based on 
the guiding-center Lagrangian, we further remove the bounce phase and obtain the bounce-averaged guiding-center 
(or bounce-center) motion. The bounce-center Lagrangian for nonrelativistic particles has been derived by Littlejohn 
8J, who at the same time constructs the second invariant Jb and the first-order correction to the second adiabatic 

10| | for relativistic particles. After 



invariant. The present work generalizes results of Littlejohn [8| and Brizard S , 
we obtain the bounce-center Lagrangian, we change the time-scale ordering of the background fields so that the 
perturbation analysis can be applied to the drift motion. We assume that the background fields vary on a time scale 
much slower than the drift time period when we derive the drift-center motion. By drift averaging the bounce-center 



3 

Lagrangian and removing the drift phase, we obtain the drift invariant Jd and the first-order correction to the third 
adiabatic invariant. 

The remainder of the paper is organized as follows. In Sec. [Til we derive the guiding-center equations of motion and 
the guiding-center Lagrangian for relativistic particles moving in slowly- varying electromagnetic fields. This Section 

n n 

generalizes previous work [J] for nonrelativistic particles and earlier work by Brizard and Chan [7[ for relativistic 
particles moving in static magnetic fields only. In addition, by introducing effective covariant potentials, we also 
simplify the relativistic guiding-center equations of motion of Grebogi and Littlejohn [6|. In Sec. IIIH we extend 
the work in Sec. Hi] and use the Lie-transform method to obtain the bounce-center Lagrangian. The first-order 
correction to the second adiabatic invariant is automatically obtained in this process. In Sec. IIV1 we assume that the 
electromagnetic fields vary on a time scale much slower than drift period and use a third Lie transform to remove 
the drift-phase dependence of the system and obtain the first-order correction to the third adiabatic invariant. A 
summary and comments on further work are given in Sec. [Vj 



II. RELATIVISTIC GUIDING-CENTER DYNAMICS 

This Section presents the guiding-center equations of motion for a relativistic particle moving in slowly-varying 
background electromagnetic fields derived by the Lie-transform method. As the first step of the Lie transform, we 
show the ordering of the background fields, then we obtain the guiding-center Lagrangian which is later used to derive 
the guiding-center equations of motion and also to obtain the bounce-center Lagrangian in Sec. IIIII 



A. Background-field orderings 



Following the work of Littlejohn 4|, we use the small parameter e = po/Lo <C 1 to order the background fields, 
where po is the typical gyroradius and Lq is the scale length of background fields. In dimensional units, e scales as 
m/q. We introduce the small parameter e by denoting the physical electric field by E p ^, and we assume that the 
Ep/j x B drift is of 0(e) compared to the particle's thermal speed % |6|, and that the background fields E p h and B 
vary on a time scale comparable to the drift period; i.e., d/dt ~ 0(e). To indicate the order of a term explicitly by 
its e factor, we set ~E p h = eE, <& p /, = e<i> and t\ — ei, where $ p h is the physical electric potential. Thus ExB~ 0(1), 
d/dti ~ 0(1), and physical results are obtained by setting e = 1. 



B. Preliminary coordinate transformation 



With the ordering of background fields given above, the particle phase-space Lagrangian one-form in slowly- 
varying background electromagnetic fields is written in terms of extended (position, momentum; time, energy) phase- 
space coordinates z = (x, p;t, W p ) as 

"1« 



-A(x,*i)+P 



• dx — W p dt — Hp da, 



(1) 



where subscript 'p' denotes particle variables and H p = H p — W p is the extended particle Hamiltonian, with H p = 
jmc 2 + qQ(x.,ti) the Hamiltonian in regular phase space. Here, the physical dynamics takes place on the surface 



Hp = 0, the guiding-center relativistic factor is 7 = y/l + \vl mc \ 2 i an d cr is an orbit parameter. 

To show the dependence of T p on the gyrophase explicitly, we decompose the relativistic momentum p according 

to 



p = p|| b + \/2m(ioBc, 



(2) 



where p\\ Q = p • b is the component of the relativistic momentum parallel to B, /^o = \p±\ 2 /2mB will be shown to be 
the lowest-order term in the asymptotic expansion of an adiabatic invariant and c is the perpendicular unit vector. 

™— j^B-— «— — 



defined by the following relations 



a = cos a e\ — sin t/ e 2 
c = — sin 9q ei — cos 9q 62 



(3) 



where c is defined by Eq. J5}, a = b x c, and (ei,&2,eg,) forms an arbitrary right-handed unit-vector set with — b. 
Substituting Eq. ([2]) into Eq. (fTJ) yields the Lagrangian written in local momentum coordinates 



--A(x,fi) +P|| b + ^2mfi Bc 



• dx 



W ? dt 



Hp da. 



(4) 



Now 7 = y 1 + 2/io-B/mc 2 + (p|| Q /mc) 2 and T p is a function of the preliminary phase-space coordinates 
(x,p|| , Ho, 9q] t, Wp). Next, we will use a Lie transform to remove the gyrophase dependence from the particle 
Lagrangian J4}. 
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C. Guiding-center Lagrangian for a relativistic particle 

A Lie transform from the preliminary coordinates z = (x,p|| , fig, 9 Q ; t, W p ) to the guiding-center coordinates Z = 
(X,f>ii, fi, 0; t, W g ) is used to remove the gyrophase dependence of T p . Here we use subscript 'g' to refer to guiding- 
center dynamics. For brevity, the steps of the guiding-center Lie transform are omitted here, but the interested reader 
may consult Ref. [121 ] for details. The resulting guiding-center Lagrangian in extended guiding-center phase-space 
coordinates (X, py , p, 9; t, W g ) is 



--A(X,t 1 )+p,|b(X,t 1 )+0(e) 
e c 



■dX 



+ e—^d9 - W e dt - H g da, (5) 
where the extended guiding-center Hamiltonian 7i g = H g — W g , with the lowest-order regular Hamiltonian 

H g = jmc 2 + q$(X,h) (6) 



= mc z ^Jl + 2^B/mc 2 + {p\\/mc) 2 + q$(X,ti). 

Here, the guiding-center coordinates are related to the preliminary coordinates and are given to lowest order in e by 
P\\ = P\\ , M = Mo, 9 = ^o, W g = W p , and 

X = x-ep, (7) 

where 



PM = - q f-^» (8) 

is the gyroradius vector in guiding-center coordinates. Note that because of the slow-time dependence, the differences 
between the guiding-center Lagrangian ([5]) and that of Ref. [7j are the electric potential and the time-changing 
variables, which give us extra second-order terms in the guiding-center equations of motion. 

D. Guiding-center equations of motion 

Having found the relativistic guid ing-center Lagrangian ([S]) , we now solve for guiding-center equations of motion 
using Euler-Lagrange equations [13|. For a Lagrangian C g , which is related to T g in Eq. §5§ by T g = C g da, written 
in extended guiding-center phase-space coordinates Z v ', the Euler-Lagrange equation is 

d tBC\ 9A =Q) (g) 



da \dZ" J dZ v 



where Z v = dZ/da. The equations of motion for t and W g are 



which indicates that t and a can be identified, and the time rate change of energy 



w g = q °—- q -K.?— , (11) 




where we replaced a by t because of Eq. (fT0| . Here, the effective potentials <£>* and A* are defined as 

<£>* i / $ 

| , (12) 

where the second term on the right side represents the covariant parallel two-flat decomposition of the relativistic 
guiding-center four- velocity [5|. 

Applying the Euler-Lagrange equation (j9]) to other guiding-center phase-space coordinates (X,p|| , p, 9), we first 
have b-X = pu/^jm), showing the parallel motion of the guiding center; secondly, 9 = e qB/^jmc), showing the 
fast gyromotion, and /i = 0, which proves that p is an invariant of the guiding-center motion (here, a dot means a 
total derivative with respect to t). Finally, the relativistic guiding-center equations for X and p\\ are 

X = ^51 + eE *x* (13) 
jmB* B* v ; 

PI, = <?E*.|! (14) 

where the effective fields (E*,B*) are defined in terms of the potentials (fT2|) as 

cr>\\ ^ 

B* = Vx A* = B + e— Vx b, (15) 

q 

and 

E* ^ V$* = E - I L ^ + mc^) , (16) 

c Oti q \ oti I 

where 

B\ = B* -b = B + e(cp\\/q)h- Vx b, (17) 

and V7 = (p/jmc 2 ) VB. Equation (fT3")) shows that the guiding-center velocity consists of the parallel motion along a 
field line, the E x B, gradient- B and curvature drifts. The curvature drift here is hidden in the first term on the right 
side of Eq. (fT3")) and the gradient drift and the E x B drift are contained in the second term. Equation (fT?)) represents 



the parallel force along a field line, which according to Eq. (fT"6)) consists of two parts: one from the parallel electric 
field and the other from the magnetic mirror force. Note that the first-order term in Eq. (|16|) gives second-order terms 
in the guiding-center equation of motion, which are important to the conservation properties of the guiding-center 
motion. 

We immediately note the simplicity of the relativistic guiding-center equations of motion (|13p and (|14p . expressed 
in terms of the covariant effective potentials (j!2p . compared to the relativistic guiding-center equations of motion of 
Grebogi and Littlejohn (GL) [6(, who used the scalar potential <I> instead of the covariant potential <&*. We recover the 
GL relativistic guiding-center equations of motion by substituting qE* = qE^ h — e(fi/j)VB in Eqs. (|13|) - (fT4|) . We 
also point out that, in contrast to Boghosian's manifestly-covariant formulation for relativistic guiding-center motion 
5|, our "1 + 3" semi-covariant formulation treats time separately from the other phase-space coordinates and uses an 
energy-like Hamiltonian (instead of the Lorentz-invariant covariant Hamiltonian) . 

If the fields are static, then Eq. (fTTj) shows conservation of energy automatically. Also, the relativistic guiding-center 
Eqs. (fT"3"l) - (fT4]) satisfy the Liouville theorem 

dBT, , .\ d 



dt V II ) d V \\ 

which ensures that guiding-center phase-space volume is conserved by the guiding-center dynamics. We prove Eq. (fTS)) 
explicitly as follows. First, we easily obtain from Eqs. (|13t - (fT4"f 

dB \\ ~ dB* db 

V- (B^X\ = c(b.VxE*-E*.Vxb) 

+ ^B*-V( 7 - 1 ), (20) 
m 

Next, we insert 



and 



«: . -Avh-'l - if, (22) 
dpu mq q dt 

^ = -Vxb, (23) 
op\\ q 



-cVxE* (24) 



at 

in Eqs. (fT^l) - ([2T|) to find that Eq. (JTSJ) is satisfied exactly, where we have set e = 1 in Eqs. (|18p 
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III. HAMILTONIAN THEORY OF BOUNCE-CENTER DYNAMICS 

To obtain the bounce-center Lagrangian, we preform a Lie transform on the relativistic guiding-center Lagrangian ([5]) 
to remove the bounce-phase dependence. This Lie transform leads to construction of the second adiabatic invariant 
and gives the first-order correction to the second adiabatic invariant directly. The nonrelativistic bounce-center 
Lagrangian has been derived by Littlejohn 8j, and the present work generalizes previous results to the relativistic 



A. Preliminary coordinate transformation 

We first drop the term e(mc/q)fid9 in the extended guiding-center Lagrangian (O, which means we are now consid- 
ering a six-dimensional system parametrized by constant-/^ surfaces. Also we separate the extended Hamiltonian TL g 
from the symplectic part of the extended phase-space Lagrangian (J5]) (as in Eq. (|26|) below) . We then perform a coor- 
dinate transformation from X to (a, (3, s), where (a, (3) are the usual Euler potentials such that B = Va x V/3 = £>b, 
and s is the position along a field line labeled by (a, (3), with b = <9X/<9s. We choose the vector potential A = aV/3, 
write 

<9X (9X <9X 

dX = —da + —d(3 + bds + —dti, (25) 
oa op oti 

and we write the (symplectic part of the) guiding-center Lagrangian ([5]) order by order as 

oo 

r g = -E e"r g „, (26) 

' n=0 

where 

L g o = ^ad/3-K s dt u (27) 
and the modified guiding-center energy coordinate 

K ^ W ^ + l a W 1 - (28) 

Eq. (|28[) introduces a change to the extended Hamiltonian 7i g = H g — Kg, where the lowest-order ordinary Hamiltonian 
H g is H g o = <7<3>* + (q/c)a dfi/dt\. It is also useful to follow Littlejohn [8( by using a 2- vector y with y\ = a, yi = (3, 
together with the two-dimensional Levi-Civita symbol r\ a b, where a, b runs overs 1 and 2. The components of rj a i, are 
given by r] n = V22 = and r) X i = -V21 = 1 - 



Before considering the first-order term T g i in Eq. (126)) written in coordinates (y, s), we make the usual assumption 



about the lowest-order motion that, with coordinates (y,ii) frozen, the bounce motion in (s,p||) space is periodic 



Thus, using the Hamilton- Jacobi theory 
to the periodic bounce motion. Then 



9]. 



13] ) we construct the action-angle canonical variables (Jbo, "0bo) corresponding 



J h0 (a,(3,fi;t,K g ) = — f p\\ds, 



(29) 



and Wbo is the lowest-order angular bounce frequency, defined by 07 bo 1 = d,Jbo/dK s . The bounce-phase angle ipbo is 
canonically conjugate to Jbo- Also, the following relation holds for the true motion: 

ds dp\\ ds dp\\ 



= 1, 



(30) 



<9i/>bo 9J b0 dJ b0 dtp b0 

since the transformation from (py, s) to (JbOiV'bo) is canonical. At lowest order, Jbo is an invariant of motion. When 
higher-order terms are included and (y, t\) are allowed to evolve, we will show that Jbo = 0(e). The symmetry of 
the unperturbed motion has been pointed out by Littlejohn 4|, and we will directly use this result to simplify the 
expression of the first-order correction to the second adiabatic invariant. 

Using the coordinates = (y, Jbo, ifibo, Wbo), the nrs * or der guiding-center Lagrangian in Eq. (|2l)|) has the 
components 



r giM = P\\ h 'Q^ = P\\ h , 



(31) 



and we will omit subscripts of Jbo an d V'bo when they themselves are subscripts. Note that in covariant form, 
b = Vs + b a \Jy a . These expressions will be further simplified with the second coordinate transformation from (py, s) 
to (J b0 , r/>bo)< 

B. Coordinate Transformation from (p\\,s) to (JbOji/'bo) 



To simplify components of r gl defined in Eq. (|3ip . we perform a coordinate transformation from (py, s) to (Jbo 5 ^bo) 
by adding a gauge term dGt, to r g i (i.e., T' gl = T g i +dGt) to eliminate the J-component (T glJ — 0). Thus we choose 
Gb to be 

JbO 



dJ' T^j^t^ipbo.J'.v)- 



Also, the ^-component = T gl ^ + dGb/dipbo becomes 



V - 



JbO 



dJ' 



9 (p\\ bj) d (pj bj) 



dJ' 



dip. 



bO 



(32) 



(33) 
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where we have used Eq. ([30)) to get the second equality and used F g u = at J b o = (since p« = 0). 
Similarly, the other components of T' gl are 



r'„ 



gla 



Jbo 



df 



and 



r; itl = / a/' 



g (P|| &o) _ d (p|| bj) 
9J' dy a 



9 (p|| b t ) 9 bj) 



0J' 



0fi 



F tl - 



Combining Eqs. ()33[) - (|35p . the new first-order guiding-center Lagrangian is 

T' gl = F a dy a + Jbo dipbo + F tl dbi. 

The Euler-Lagrange equation [see Eq. ©] for J b o obtained from T' g = F g o + eT' gl = C' g da is 

rfJbo dF a . dF h . 2 , 
at Olpbo Olpbo 



where to lowest order 



dH m 



Va = e-Vab- 



•fau <~x 

q oyb 

For later use, we now write the first two terms of the guiding-center Lagrangian [Eq. 
(a, j3, Jbo, ipbo', ti, Kg) with the prime dropped, 

r g o = —adf3 — K g dti, 

r g i = F a dy a + Jbodipbo + F tl dti, 



and the lowest order Hamiltonian is 



7~t g - -ffgo - K t 



(34) 
(35) 
(36) 
(37) 
(38) 

in coordinates 

(39) 
(40) 

(41) 



With these coordinate transformations and Lagrangian, we do a Lie transform to remove the bounce-phase depen- 
dence from T g and obtain the bounce-center Lagrangian r b . 

C. Lie transform in extended phase-space coordinates 

The bounce-center dynamics are obtained using the Lie transform in extended phase-space coordinates. The 
bounce-center Lagrangian and Hamiltonian are constructed order by order 



Tb = r b0 + eF bl + e 2 F b2 + 
Wb = Ti-bo + e^bi + £ 2 ?ib2 



(42) 
(43) 
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where the terms on the right-hand side of Eq. (|42|) are 

r b0 = r g0 , (44) 

r b i = T g i - ii ■ fl g0 + dSi, (45) 

Tb2 — ~*2 • ^gO — *1 ' ^gl 

+ t -±.d(i 1 -Sl g0 )+dS 2 , (46) 

and the first two terms in Eq. (|4"3")l are 

Wbo = W g o, (47) 

Ti-bi — W g i — <?i • dHgo- (48) 

The term i„ • Q g = g£(Q e ) ^dZ" in Eqs. 05} - (gBJ) and the term g„ • dft g = g£ ■ dH g /dZ^ in Eq. (gSJ) are expressed 
in terms of the n t,l -order Lie-transform generating vector g n and gauge function S n , where 

(^W-[^^] = ^-|^ (49) 

is the Lagrange bracket between Z* 1 and 

D. Bounce-center motion in coordinates (Y , Jb,ipb',t, Kb) 

Following the Lie-transform procedure described in Eqs. l|4"4"] ) - (|38"]l . we first have the lowest-order Lagrangian and 
Hamiltonian 

T b0 = ?-Y a dYp - K b dti, (50) 
c 

n b0 = H bQ -K bl (51) 

where (Y a ,Yp) represent the bounce-center coordinates Y, and H b Q has the same functional dependence on the 
bounce-center coordinates (Y, J b ,t) as H g o on the guiding-center coordinates (y, J b o,t). 
The first-order bounce-center Lagrangian (|45[) then becomes 

T b i = {-gfaabo + F a )dy a + J b0 d^bo + [F tl + g?)dh, (52) 

where Q a bQ = — (<z/c)f? a b and we choose S± — in Eq. (|4"5"|) . Requiring Tbi a = and T b i tl — gives us 

gl = --VabF b , (53) 

q 

,K 



9i = -Ft,- (54) 
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The first-order Hamiltonian then is given by 

n bl = -g a 1 ^-g J 1 u b0 + g^, (55) 



since H g i = in Eq. (|48[). With Eqs. ([53]) . (|54|) . and the requirement Tib = 0, where a tilde in this section denotes 
the bounce-phase oscillatory part, we have 

7i = — (-F tl + -Va b F b ^) , (56) 

and 

Wbi = -Vab(F b )^ - ( 9 ()uj b0 - (F tl ), (57) 

n 

where (• • • ) denotes a bounce-phase average. It has been shown in Ref. IS] that F a and F tl are odd in ipbo, and thus 
we have (F a ) = and (F tl ) = 0. Equation (j57|) then becomes 

Wbi = -<<?i>bo- (58) 

To obtain gf and the bounce-phase averaged part of g{ needed in Eq. (|58|) , we need to go to the second-order Lie 

transform of the Lagrangian. The V'bo part and the Jbo part of the second order Lagrangian r b 2 are 

dS 2 j 1 c dF a 

r b 2V = 9i - ^-VabFb-^—, (59) 

dip b o 2g di/>bo 

r b 2./ = +g\ [- »-VabF b ^-— . 60) 
d Jbo 2 g d Jbo 

To make Jb the exact invariant, we require that T^ip = 0. Taking the bounce-averaged part of Eq. (f59|) and using 
(S2) = 0, we have (g() = 0, since dF a /dipbo is even in i/j^o, and thus (Ft dF a /d4>bo) = 0. This result indicates that 
Eq. (|58l) becomes H b i = 0. 

The bounce-phase dependent part of Eq. (|59[) is solved to give the gauge function 



S 2 = / dipbo 9i + -z- r labF b7 — 
J \ 2q dtpboj 



(61) 



where g{ — g( is given in Eq. (|56|). 

Inserting S 2 into Eq. (|60p and requiring L b 2j = 0, such that Jb and ipb are exact conjugate coordinates, gives us 



v, 1 c „ dF Q ds 2 



2 g d J b0 d Jbo 



Thus we obtain the bounce-center Lagrangian and Hamiltonian 



T b = llY a dYp + J b di> b - K b dt, (63) 
e c 

H b = H b0 + O(e 2 ), (64) 
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where the bounce-center coordinates (Y, Jb, ipb] t, Kb) are given by 

Y a = y a - e-r) ab F b + £>(e 2 ), 

Jb = Jbo + esf + 0{e 2 ), 

A = ipbo + egf + 0(e 2 ), 

K b = K g -eF tl +0(e 2 ), 



(65) 
(66) 
(67) 
(68) 



with time t an invariant under the transformation. Note from Eq. (165[) that F a in bounce-center dynamics plays a 
role similar to the gyroradius vector p in guiding-center dynamics [Eq. |7|)]; i.e., F a may be interpreted as a "bounce 
radius" 2-vector. Also from Eq. (jB"5)) , i 7 ^ is the oscillatory part of the guiding-center energy coordinate K g . 
The bounce-center equations are then 

dUb 



Ipb = 
Jb = 

K b = 



q la " dY b > 
dHb 
dJ b ' 

0, 

dHb 



(69) 
(70) 
(71) 
(72) 



dt ■ 

Thus we see that Jb is the exact invariant for the bounce motion. 

In Eq. (|56|) . <?/ = g( denotes the first-order correction to the second adiabatic invariant Jt,o- This first-order 
correction can also be directly obtained from Ref. [lj, where Northrop shows that 



dJ\ 



dt 



b0 -1 
W b0 



l((a)P-a(0)) + (K g -(K e ) 



(73) 



written using our notation. Since 



dJb _ rfJbo ^dJbi 
dt dt dt 



= 0, 



(74) 



and to lowest order, we have (dJbo/dt) = 0, thus 



Jbi = -w. 



bo 



y bo 



r/J, 



bO 



rfV'bO 



#b0 

| ((d)j9 - a0) 



+ [K g -(K g } 



(75) 



Since we have shown that Jbi (i.e., g-) is purely oscillatory, we have Jbi = Jbi- Littlejohn 8] has shown, for the 
nonrelativistic case, that the right hand side of Eq. (fTS"]) is equal to the right hand side of Eq. (|5"6"|) , but this result 
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by Littlejohn 



also holds for the relativistic case because the equations have the same functional form. We also note, as pointed out 
8J, that the Lie-transform approach is more straightforward than the method used in Ref. 



14j . which 



derives the first-order correction to the second adiabatic invariant for nonrelativistic particles in a static magnetic 
field. 



IV. HAMILTONIAN THEORY OF DRIFT-CENTER DYNAMICS 

Starting from the bounce-center Lagrangian (|63p . we now derive the drift-averaged bounce-center Lagrangian, or 
the drift-center Lagrangian. Similar to the analysis given in Sec. IIII1 this procedure leads to the first-order correction 
to the third adiabatic invariant automatically. 

To apply the adiabatic theory to the drift motion, electromagnetic fields must vary on a time scale much slower 
than the drift period; i.e., d/dt ~ e 2 . We start from the bounce-center Lagrangian (|63[) with term Jbdipt, dropped, 
which means we are now considering a two-dimensional motion parametrized by the constants J g and Jb- We set 
t 2 = e 2 t and the resulting bounce-center Lagrangian is 

T b = -ladp - \x b dt 2 = \ (eadp - K h dt 2 ) , (76) 
e c e z e z 

where we henceforth use Y = (a, 0) and replaced qa/c by a in the last expression. 

We now make the usual assumption for the lowest-order motion that in a static field, or with parameter t 2 frozen, 

the orbit of the particle is closed and hence the drift motion of the particle is periodic l|. Thus the coordinates 

(a, (3) play a role in drift-center dynamics similar to that of as the coordinates , s) in bounce-center dynamics. The 

Hamilton- Jacobi theory again gives us the action-angle variables from canonical coordinates (a, (3) as 

Jdo(K b ,t) = i- j ad/3, (77) 

and oJjq 1 = dJdo/dKt, is the lowest-order angular frequency of the drift motion. Here we use 'd' to represent drift 
motion variables. The canonically-conjugate coordinate of Jdo is the drift phase 7/^0 • The change from coordinates 
(a, (3) to ( Jdo, V^o) is canonical, thus we have 

da d(3 da 3(3 =i 
dJdo dipdo dipdo dJdo 

which is also valid for the true motion. 

For the lowest-order motion, J^q is a constant, but with time t unfrozen and higher-order terms included in the true 
motion, J^o is no longer an invariant for the drift motion and it will be shown that dJ^o/dt = 0(e 2 ). Thus we first 
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do a coordinate transformation from (a, 0) to (Jdo, V'do) and then use a Lie transform to construct the true invariant 
Jd for the drift motion. 

A. Preliminary coordinate transformation 

Similar to the construction of the bounce-center dynamics in Sec. IIII1 we first change coordinates from (a, (3) to 
PdOi^do)- Substituting 



dp dp dp 

oJdo oipdo dt 2 



which is similar to Eq. (|25"]l . into Eq. ([76]) gives 



e T b = ea— — d J d0 + ea— — dip d0 



which gives the lowest- and first-order bounce-center Lagrangians 



Kw — ea 



dp_ 



dU 



(79) 



(80) 



Tbo = —Kbdtz, 

Tbi = a——dJ do + a——d-ip d o + a—dt 2 . 
dJ d0 oip d o dt 2 

Similar to Eq. (|32[) . we perform a gauge transformation on Tbi; i.e., r' bl = Tbi + dGd, such that 

r'bi,/ = and r bl ^ = J d0 , 

where we have again omitted the subscripts of Jdo and V'do when they themselves are subscripts. From Eq. 
choose Gd as 



(81) 
(82) 



G d 



T b udJ' + f(tpdo,t 2 ), 



(83) 



, we 



(84) 



which is similar to Eq. Q32p . and /(V'do 1*2) is a function that is determined from the condition 

SG d 



Oil), 



Since Tbiw, can also be written as 



Tbiv 



do 



dT 



rbii/i = Jdo- 



(85) 



blV 



the equation for /(V'dO;^) becomes 



dJ' 



df _ dp 
a- 



—dJ' + r b i^| Jbo=0 , 



(86) 



<9V>do dijj, 



do 



(87) 



Jbo=0 
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Note the difference between Eqs. ((84)) and (|32| . since we generally do not have rbi^|j bo=0 = ^- Finally the gauge 
transformation (|84|) yields the new ^-term 

= F t2 {J M ^ Ml t 2 ). (88) 
Now we have our zeroth- and first-order bounce-center Lagrangian 

r b0 = -K b dt 2 , (89) 
r bl = Jdo#do +F t2 dt 2 , (90) 

where we have dropped the prime, with the extended Hamiltonian 

n b0 = H hQ - K h . (91) 

From Tbo and Tbi, we obtain the Euler-Lagrange equation for J^o 

d J^ = ^ + 0{e% (92) 
at dip 

Comparing Eqs. (|89"|) - (f9"Tj) with Eqs. (f3"9")) - (|41[) . we find that the bounce-center and guiding-center equations are very 
similar, except that in Eq. (|89[) we do not have the d(3 term and in Eq. (|90[) wc do not have the dy a term. Thus 
the Lie transform from the bounce-center coordinates to the drift-center coordinates will be very similar to the Lie 
transform from the guiding-center coordinates to the bounce-center coordinates. 

B. Lie Transform from ( J i0 , tjj d0 ; t, Kb) to (Ja,tpd',t, K d ) 

The lowest-order drift-center Lagrangian and Hamiltonian are given by T^o = Tbo and Hda = Ti-bo- The first-order 
Lagrangian is given by 

rdi = Tbi - i\ ■ Obo + dS\. (93) 
Choosing Si = and substituting ij • fibo = —Q\dti and Tbi from Eq. (|9H|) into Eq. (|9"3"|) . we obtain 

Tdi = Jdo#do + (Ft 2 + 9i)dt 2 . (94) 
The first-order Hamiltonian is then given by 

Hdi = Hbi - 3i • dHbo = -g'l^do + 9i , (95) 
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where we have used Hti = 0. Note that Eqs. (|94| - (|95)) look similar to Eqs. (|52| and (|55|) and we omit the remaining 
details here. 

The first-order coordinate generators from the above Lie transform are 



9( = 1 r±=-u>£K, (96) 
dSo 



9l = ~' (9?) 



Jl 

9i = -Ft,, (98) 
where a tilde in this Section indicates the drift-phase oscillatory part. Thus the overall coordinate transformation is 

J d = Jao-eu^Kz +0{e 2 ), (99) 

= fao + egf + 0(e 2 ), (100) 

K d = K b -eF t2 +0{e 2 ). (101) 

The drift-center Lagrangian written in coordinates (Jd, ipd] t, K$) is 

r d = - J d d^ d - K d dt, (102) 
e 



and the drift-center Hamiltonian function is 



H d =^ d0 + eHdi +0(e 2 ), (103) 



where Hdo — Hto and TCdi = —((Ft 2 )), with a drift-phase average denoted as ((••■})■ 
The drift-center equations of motion are obtained from the Euler-Lagrange equations 

Jd = 0, (104) 

* " §£• <105> 

* - it- « 

Equation (|99|) shows the first-order correction to the third adiabatic invariant. We can also write the oscillatory 
part of F t2 in another form by using Eq. (|92[) and 

dJdo dJ d0 dJdo ■ . . 

^r = ^r + dK- b Kb > (107) 

which gives that 



f,, = ; #, 



/"^(tt + s;*)' (108) 
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where we have set e = 1. Northrop [l| has shown that the first term on the right side of Eq. (|107[) can also be written 



With d Jdo / dKt, = Lo d ^ , we find 



n 



a j. 



do 



dt 



F t2 = c^ 1 / dVdo [Kb ~ «^b» 



(109) 



(110) 



Thus, Ref. [l[ implicitly contains the first-order correction term to the third adiabatic invariant (see Eq. (3.80) on 
page 64 of Ref. |lj), but Eq. ([99]) is an explicit expression. 



V. SUMMARY AND DISCUSSION 



In this work, we have presented the Hamiltonian theory of adiabatic motion of a relativistic charged particle and the 
derivation of the first-order corrections to the second and third adiabatic invariants. The background electromagnetic 
fields vary on the drift time scale when we consider the guiding-center motion and the bounce-center motion. The 
effect of these time- varying background fields on the guiding-center motion is shown by the extra terms in the guiding- 
center equations ()13|) and (|14j) . compared to the guiding-center equations in Ref. 7[. The first-order correction to 
the second adiabatic invariant of a relativistic particle is then shown in Eq. (|66|) . To apply the adiabatic analysis to 
the drift motion, we assume that the background fields vary on a time scale much smaller than the drift period. The 
first-order correction to the third adiabatic invariant is shown in Eq. (|99p . 

This work simplifies previous work on relativistic guiding-center motion, generalizes previous work on bounce-center 
motion for a relativistic particle in time-varying fields, and extends previous work on drift-center motion using Lie- 
transform perturbation methods in extended phase space. These results are especially useful in space plasma physics, 
where adiabatic theory is the foundation for modeling and understanding the dynamics of magnetically-trapped 
energetic particles. 

The hierarchy of the adiabatic motions in this work may be shown as follows 



(x,p;t,W p ) { 



(X.piijt.Wg 



d (*,K"d) 
(a,(3;t,K b ) -5-> ! 
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where J g is related to the first invariant /i by J g — (mq/c) /j, and ip g = 8. The first arrow (g) thus indicates the 
gyro-phase average process, the second arrow (b) the bounce-phase average and the third arrow (d) the drift-phase 
average. 

In this paper we have shown how first-order corrections to adiabatic invariants can be obtained using Lie-transform 
methods. Alternatively, the oscillatory part of the first-order correction to an adiabatic invariant can be obtained as 
follows. Diffcrenting the exact invariant 



Jk = Jko + eJki + ■ ■ ■ , (HI) 



where k can be g, b or d, to lowest order gives 



dt dt dipko 

Since Jko satisfies the necessary condition 

/ dJko \ 



dJ k d,J k0 , d,J k i Hm\ 
ewko-r- 1 = 0. (112) 



= 0, (113) 

k 



\ dt j 

where (• • • )k denotes the fast-angle average canonically conjugate to Jk, and we obtain the oscillatory part of Jki as 



eJki = -w ko x J -^T*Ako- (114) 

The phase-independent part of Jki can be obtained by using the Lie-transform method. 

The use of Hamiltonian theory in describing adiabatic motions results in equations that satisfy energy conservation 
for time- independent fields and preserve phase-space volume naturally, in contrast to the results of Northrop 1] . These 
conservation laws are very useful for checking numerical accuracy in simulations. Based on this work, fluctuations of 
electromagnetic fields can be added to the background fields and equations of motion in the presence of electromagnetic 



waves can be derived, as in Refs. [7|, [9(, and 
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